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Abstract 

The modification of the width of rho mesons due to in-medium decays and collisions 
is evaluated. The decay width is calculated from the imaginary part of the one-loop self- 
energy at finite temperature. The collision width is related to the cross sections of the p + 
pion and the p + nucleon reactions. A calculation based on an effective Lagrangian shows 
the importance of including the direct pn — > pn scattering which is dominated by the ai 
exchange. A large broadening of the spectral function is found, accompanied by a strength 
suppression at the pole. 



PACS : 25.75.Dw, 14.40Cs, 21.65.+f, 13.75.Lb 



I. INTRODUCTION 



It is generally believed that ultrarelativistic heavy ion collisions produce matter far from the 
ground state and provide conditions of high temperature and/or high density suitable for investi- 
gating modified particle properties in the medium. The vector mesons, e.g. the rho, are particularly 
important components of the matter produced in ultrarelativistic heavy ion collisions. It has been 
speculated that the restoration of chiral symmetry induces a mass shift of the p meson fij. The 
direct dilepton decay channel, p — > e + e~ , appears to be an ideal probe for studying the in-medium 
properties of the p: a possible mass shift, a modified dispersion relation, or a width modification 
Unfortunately, due to the interplay between those effects, it seems to be difficult to iso- 
late the different effects from invariant mass measurement alone. The in-medium mass shift and 
broadening of the width of a particle have been related to the real part of the forward scattering 
amplitude on the constituents of the medium and thus to the corresponding cross sections. These 
results indicate a small mass shift but a large increase of the width both in intermediate and high 
energy heavy ion collisions Pj9|,|l"c|1 . The goal of the present paper is to verify these results using 
slightly different techniques. We will point out the differences as we proceed. 

Recently, dilepton measurements at the CERN/SPS energy have received quite some attention 
The e + e~ mass spectra (in S + Au collision at E/A = 200 GeV and in Pb + Au collision 
at E/A = 160 GeV preliminary data) seem to indicate a /3-peak suppression in central collisions. 
The absence of p-peak may be due to a large width of the p meson as compared with its free value. 
Of course, a confirmation of such a conjecture must include a dynamical simulation for the space 
time evolution of the colliding system. In any case, it is important to quantitatively address this 
issue independently. 

The in-medium effects on the p meson have been widely investigated using finite temperature 
field theory at the one- loop level p.l2,13||) but the most important source for the width modification 
occurs at the two-loop level ||. The present paper focuses on the p meson width modification up 
to the two-loop level. The pole shift of the p, as emerged from many-body calculations is small 
IT( 
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The in-medium properties of the p are due to its interaction with the constituents of the 
medium. In the hadronic matter produced in ultrarelativistic heavy ion collisions, pions and 
nucleons are important constituents, so we take into account both pn and p N reactions in our p 
width calculations. Interactions with other hadrons are expected to increase the width above the 
values found here. 

The paper is organized as follows: In section II, we present the general formula? for the collision 
rate of in-medium p's. We present a calculation of the thermal averaged in-medium width. In 
section III, we discuss the cross sections for pn- and /?N- reactions from which the collision rates 
are determined. The width modifications due to pix- and pN- reactions are calculated at finite 
temperature and density in section IV. Section V concludes with a summary. 



II. FORMALISM 

The principal elementary reactions which tend to thermalize p's in hot and dense matter are 
the channels p\ 1x2 - > P3 ^"4 and/or p\ N2 — ► p 3 N4. A p with arbitrary momentum pi (and energy 
uS) can be captured by the thermal background. The rate for this direct capture |l4| is: 

r d (u) = ^f dnn 2 (E 2 ) (1 ± n 3 {E 3 )) (1 ± n 4 (# 4 )) |A1(12 - 34) | 2 , (2.1) 

A p with momentum p\ can also be produced from the thermal background by the inverse 
reaction: p 3 tta — ► p\ 712- This inverse rate (omitted in many similar treatments) is 

T ^) = ^j dQ C 1 ± ME2)) n 3 (E 3 ) ni {E 4 ) |A4(34 ^ 12)| 2 , (2.2) 



where 



and 



dn = dp 2 dp 3 dp 4 (27r) 4 S(pi +P2-P3- Pi) , (2.3) 



d? Pl 

dPi 



(2vr)3 2Ei ' 

ni(Ei) are Bose-Einstein or Fermi-Dirac distribution functions, 'db' corresponds to Bosons and 
Fermions, respectively. 



The total collision rate for bosons is the difference 

r coll H = r d H-r i (a;). (2.4) 

This expression can be derived from quantum field theory [15]. From finite temperature quan- 



tum field theory, the width of any particle, r(w), is given by the imaginary part of its self-energy, 
Imn(w), 

Imn(w) = -uT(u) , (2.5) 

which provides a link between field theory and statistical mechanics. 
The ratio of Td(uj) to Ti(u) is an universal function of u flT6|| , 

= exp( W /T) , (2.6) 

because the single-particle distribution functions satisfy 

1 + rasM / //ti\ /-p x 

r—- — = expuj/i ) (Boson) , 

n B {uJ) 

l-n F (w) 

r— — = expfu/i ) (.bermion) . (2.7 ) 

Thus the total collision rate is given by: 



-icoll 



(l-e-^)r,H . (2i 



In the simplest scenario of equilibrated hot and dense matter, we can calculate the thermal 
average of the collision rate: 

pooll = Jd 3 Pir ro \^) = Jd 3 Pl T d e-^ 



Using ( |2.1| ) and (2.9), the thermal rate can be written as 



r° 011 = 1 2 f ds T \(s,mj,ml)K(^,m 1 ,m 2 ,T,Li i )<T 1 2^34(s) , (2.10) 
pi J so 2(27r) 4 Vs 



where 



Pi = Mi I («") (2- 11 ) 



is the number density of the p's we considered, ai2-+34(s) is the cross section of pit (or pN) 
collisions, and sq = (mi + m2) 2 , Mi are the spin-isospin degeneracy factors. The function K 
assumes different forms for different statistics [jT^ . Below we apply classical Boltzmann statistics, 
where 

K = K 1 (^/T)exp[(fi 1 + p l2 )/T], (2.12) 

with K\ being the modified Bessel function. p\ , p 2 are the chemical potentials of p and ir (or 
nucleon) respectively. From ( p. 10 ) the collision rate of p in the hot and dense matter can be 



obtained if the cross sections of the p it- or pN- processes are known. 

III. pir AND pN CROSS SECTIONS 

For pir collisions, we use both the scattering amplitude of the elementary processes method 
and the resonance contribution approximation to calculate the cross sections a pn . We consider the 
scattering amplitudes for the elementary reactions in Fig. 1 which can be obtained by cutting the 
two- loop diagrams for the p meson self-energy as demonstrated in Ref. [18]. 



To determine the p N cross section a p N, a method based on the resonance model is used to get 
a lower estimate. 

A. pir cross section 

The standard Breit-Wigner formula for p it cross section is 

1 r (Vs - m R y + T|j/4 

Here 'R' refers to the intermediate resonance, and -Br is the branching ratio of its decay into n p. 
yfs is the total center-of-mass energy, q is the three-momentum of the p meson in the resonance 
frame, 

q = -L\i/*(s,ml,rn 2 p ), (3.2) 



with the kinematical triangle function X(x,y,z) = x 2 — 2x(y + z) + (y — z) 2 . F s and F{ are spin 



and isospin degeneracy factors: 

2I R + l 



F 



F x 



2Jr + 1 



(3.3) 



(2J ff + l)(2/ p + 1) ' ^ s (2 S?r + l)(2s p + 1 
We take into account the following resonances: 0(1020), ai(1260), o 2 (1320), u/(1420), vr(1670) 



with the properties listed in [19|. 
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FIG. 1. The Feynman diagrams, (a) for p ir — > p ir reactions and (b) for p ir — * p"ir 
reactions. 



The pion scattering amplitude on any hadronic target vanishes at zero pion energy in the target 
rest frame in the zero pion mass limit (Adler's theorem). Following Ref. |J one can include an 
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additional factor to the resonance contributions to a p7T , 

(3.4) 



2 2 \ 2 

s — m p — \ 



I 2 2 2 I ' 

\ m R — — m£ y 

which is normalized to 1 at s = m R ; for s > m R this factor is taken to be one. The resulting 
resonant cross section a pn is shown in Fig. 2(a) (dashed curve). 

The cross section for pir scattering can be calculated by using an effective Lagrangian for the 
hadronic interaction pix — * pir. The scattering proceeds through s- and t- channels respectively, 
where the it, p, u>, (f>, a\ and u/(1420) might be intermediate states. Those processes are constructed 
from Trirp, irpuj and ixpa\ vertices. The effective interaction Lagrangians we use are the following 



Cpwn = igpmr(n x d^ir) ■ pV , (3.5a) 
£p P TTn = 9%^ x Ptx) ■ (if x pt") , (3.5b) 
£p*u> = g f mu>e >u " rr d ll u l/ d ( rp T ■ tt . (3.5c) 



The Lagrangian £ p7ra i is too long to be presented here [21,22] but for the physical a\{k) — > 
p(q)ir(p) decay, for example, the vertex function is 

IV = i{fA9piv + gAQ^K + h A q^q v ) , (3.6) 

where 

fA = —/= ~ViQ 2 + (Vi ~ V2)k ■ q , 

9a = -^(vi -m), 

h A = —J?! . (3.7) 

The third term in ( |3.6| ) actually does not contribute to the decay width. We chose parameters 
which reproduce the measured a\ — ► pir decay width and D/S ratio (g = 6.45, r/i = 2.39 GeV -1 , 
r\2 = 1-94 GeV -1 ). The coupling constant g pn7T = 6.05 is determined from a fit to the free p — > irir 
decay width, g pn0J = 12 AO GeV -1 is determined by the uj — > pn — ► ir^y decay [ 20 1 . The coupling 



constant g piTU} i = 3.55 GeV" 1 is determined by the uj' — > pn decay. Note that we believe that it is 
important for our effective Lagrangians to generate good phenomenology at energy scales relevant 
for the applications considered here. 

The Feynman diagrams for pir scattering are given in Fig. 1. The cross section is then calculated 



as 



dt\M\ 2 , (3. 



167rA(s, ml,m 2 p ) j t _ 
where 

t± = ml + vn 2 - ^-[(s + ml- m 2 p ) (s - m 2 + m 2 p )} ± ^A(s, m 2 , m 2 p ) . (3.9) 

For p° 7r^~ — > p° tt^~ scattering, the squared amplitude is obtained from the processes shown in 
Fig. 1(a) with the appropriate interference contribution: 

1 



\ M \ = (o xivo ..A Mi+Mz + Mz + Mi + Mtf, (3.10) 
(25^ + l)(2s p + 1) 

which is sum over the spins of the initial states as well as over those of the final states. The 
p° -k^ — ► p^~ 7r" scattering is also considered, The squared amplitude has the same form as ( p.10 ) 



but no contributions from A^2 and For p°7r° scattering, Fig. 1(b), \M\ can be calculated 

similarly. 

Due to the substructure of hadrons, for vertices in the t-channel Feynman diagrams we use 
form factors 



A 2 



,2 



Fa = -jr^f , (3.11) 

where a indicates the species of the exchange particle. A = 1.8 GeV is taken for all vertices. 
The numerical integration for the reaction pn — > tt — > pir (t-channel) results in a singularity. We 
regulate it with the effective approach of Peierls . This technique offers a pragmatic solution to 
the possible unitarity violation at high energies that can result from the use of tree-level diagrams. 
A more complete solution would perhaps comprise a K-matrix multichannel calculation. However, 
such an endeavor at finite temperature and finite density also carries ambiguities of its own. 
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After computing the amplitude, the tt p scattering cross section is obtained from ( |3.8|) . The 
result is shown in Fig. 2(a) (solid curve). The resulting pn cross sections indicate that the difference 
between the resonance model and the scattering amplitude method becomes large at high energies. 
The resonance model may be valid only at low energies. Our effective Lagrangian technique also 
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FIG. 2. The cross sections, (a) for pn — > pn reactions, and (b) for p N — ► p N reactions. 



B. pN cross section 



The p N cross section can be calculated within the resonance model as M 



apN = ^ (2jR + 1) _0_ . (3 . 12) 



Here 



q =-^=X l /\ s , m %,ml), (3.13) 

and y/s Q = m p + mjv- The summation is performed over all baryonic resonances with masses 
above the pH threshold and below 2200 MeV as quoted in The result of the pN cross section 



calculated from ( 3.12 ) is shown in Fig. 2(b) (dashed curve). 

The pN cross section obtained within the resonance model may be valid only for low energies, 
while at high energies one should calculate it from the quark model. Furthermore, it has been 
shown recently that the resonances below the p N threshold, like the N*(1520) play an important 



role in the p-nucleon dynamics [24|, so in the baryonic sector we adopt the pN cross section as 
determined by Kondratyuk et al. Ref p5[ , which includes also resonances below the pN threshold. 
The numerical result is shown in Fig. 2(b) (solid curve). 
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IV. BROADENING THE WIDTH OF THE RHO AT FINITE TEMPERATURES 

AND DENSITIES 



The modified properties of hadrons in hadronic matter, e.g. the width increase of the particles, 
are due to their interactions with the constituents of the medium. We assume that ultrarelativistic 
heavy ion collisions generate a hot and dense hadron gas of pions and nucleons and consider the 
p mesons interacting with this thermal gas. We realize that this is a simplification of the complex 
dynamics which will be generated in all collision of heavy ions. Our results should be taken as 
lower estimates of the realistic effects. 

As discussed in section II, the rho meson width is related to the imaginary part of the rho 
self-energy. The total width of the p meson is the physical decay width plus the collision width, 

ptotal = pdecay + j^coll. _ ^ ^ 

The decay width is related to the imaginary part of the rho self-energy at the one-loop level. 
The collision width is related to the imaginary part of the rho self-energy at the two-loop level. 
Let us first calculate the decay width of the p mesons at finite temperature. 

A. The decay width of the rho at finite temperature 

To evaluate the p meson self-energy from the tttt loop we consider the pnir interaction La- 
grangian from the previous section. Using the Feynman rules of thermal field dynamics, the 
temperature dependent rho meson self-energy is [jl^] 

IW?, T) = J d 4 k {-2g^ v 5(k 2 - ml)n^{k) + {q + 2k) + 2k) v - 



5(k 2 -ml) n . 5((k + q) 2 -ml) . 

rn^k) + >2 v, n*{k + q) 



(k + q) 2 — m 2 71 k 2 



m 



(4.2) 



Because of the lack of Lorentz invariance in the medium p,12,26|], the p self-energy at finite 
temperature separates into a transverse part, TlT(q,T), and a longitudinal part, Hi(q,T). In the 
rho meson rest frame, 

n(w,g^0,T) =n L (w,g^0,T) =U T {uj,q^ 0,T) (4.3) 
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thus the real part of the self-energy is 
ReU(uj,q -> 0,T) = 



4k 2 



n n (E k 



(4.4) 



where E k = \Jk 2 + m^. The real self-energy will lead to a mass shift of the p at finite temperature, 
which is determined by the solution of 



u? -m 2 + ReU(u;,q -» 0,T) = 0. 



It is known from previous calculations that this mass shift is small 



(4.5) 



The imaginary part of the self-energy can be obtained by using cutting rules |16|] , 



lmU(co, T) = ^ / dk^- (2n^(E k ) + 1) [S(u + 2£ fc ) - 6(u - 2E k )] . 

67T J E k 

From ImIl(u;,T), the corresponding width of the p — > 7T7T decay follows to be 



(4.6) 



r^ eca y(T) 



Q 2 
48vrw 2 



(u, 2 -4m 2 ) 3 / 2 



(4.7) 



This decay width varies with the temperature as shown in Fig. 3. Obviously the decay width 



increase^, 



/'tI 



with temperature as well as with the pion chgmreal potential 



w ■ o.-^r,, = r';, ,,,a > i/'i rf ,a > (o) = 25 




;V, while lor =p 
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135 MeV, Ar p 70 MeY. 




-T (GeV) 



-T (GeV) 



FIG. 3. The p — ► 7T7T decay width as 
a function of temperature T at two fixed 
values of the pion chemical potential. 



FIG. 4. The p collision rate vs. the tem- 
perature due to p 7r interaction in two dif- 
ferent approaches. 



B. Collision rate of rho at finite temperature and density 



From (230) and with the cross sections computed in the previous section, the collision rates of 



the p mesons can be calculated numerically. 

First the rho collision rate is calculated due to pir reactions in the finite temperature meson 
medium. The contributions from p ^ — > p ^ , p ^ — ► p^ir and p°7r° — ► p°7r° processes are 
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taken into account. Fig. 4 shows the the collision rates as calculated from the scattering amplitude 
calculation and from the Breit-Wigner resonance approximation at two values of the pion chemical 
potential, p,^ = and p n = 135 MeV. For the same temperature and chemical potentials and using 
the pn cross section from the scattering amplitude method, one can get a larger collision rate of p 
mesons. 

The collision rate from p N reactions is plotted in Fig. 5. In Fig. 5(a), we use T = 100, 150 
MeV and show the rate as a function of the nucleon density ratio. In Fig. 5(b), we show the rate 
as a function of the temperature for two values of the nucleon chemical potential, pn = 0.4, 0.6 
GeV. From Fig. 5(a), we can see the collision rate increasing with the nucleon density for a given 
temperature, but decreasing with temperature for a fixed nucleon density. The latter decrease is 
due to the fact that the nucleon chemical potential decreases with temperature at a fixed nucleon 
density. Again we find that for given temperature and chemical potential the pure resonance 
contribution yields a lower collision rate than the resonance + quark model. 
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FIG. 5. The collision rate of the p due to p N interactions, vs. (a) the nucleon density for vary 
temperature, (b) the temperature for vary baryon chemical potential. 



We employ the resonance model modified to include the resonances below the pN threshold, 
such as N*(1520), as described previously. The total pN cross section can be evaluated as a function 
of the invariant collision energy i/i and the invariant mass of the p meson M p . In this case, the 
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collision width has to be integrated over the spectral function of the p meson as 



ds- 



*fs-m N 



dM p M p A(M p ) \(s, Mlm 2 N )Ka pN (s, M p ) , (4.8) 



Pl Js > 2(2^) 4 ^ Am. 
here s' = (2m n + m^) 2 . A(M p ) is the spectral function of the p meson in free space taken as 



A{M p 



m p T p {M p ) 



1 



7T (M2 - m 2 p ) 2 + m 2 p Tj 



(4.9) 



where m p = 770 MeV, Y p {M p ) is the mass-dependent width of the p meson [25]. For the same 
process of pN scattering, the collision width of the p meson in this calculation is about 10% larger 
than that from the on mass-shell p meson calculations. 

In summary, our results indicate that at high temperature and low nucleon density, the collision 
rate is dominated by p + meson reactions, while at low temperature and high nucleon density it 



is dominated by p + nucleon reactions. For example, when T = 100 MeV, p n = 0, P n 



Po, 



= coll. 



= coll. 



12 MeV, while F pN 



180MeV. 



For conditions typical of high energy heavy-ion collisions |,||, T = 150 MeV, p n = and 
/ijv = 0.4 GeV, (pn = 0.39/9o ) the total width of rho mesons is 

ptotal _ pdccay _|_ p co11 - _|_ p 

« 176 + 58 + 63 = 297 MeV . 
This dramatic increase of the p width is most certainly important. 
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FIG. 6. The spectral function of the p vs. the invariant mass M of the p for fixed values of 
momentum as indicated, (a) In pure pion gas. (b) In the matter of pions and nucleons. 

The spectral function of rho meson is given by 

2Imn 

P( ' ~ ~ (u 2 - q 2 - m 2 - Ren) 2 + (Imll) 2 ' ( ' 10) 

where q is the three-momentum of the p meson. The vacuum part of the self-energy II, which is 
given in |2g|, can only depend on the invariant mass M = y/oj' 2 -q 2 . We plot the p meson spectral 



function in fig. 6(a) for a pure pion gas at a temperature T = 150 MeV and in fig. 6(b) for the 
matter of pions and nucleons at T = 100 MeV and nucleon density p^ = po = 0.16 fm~ 3 . We find 
a broadening of the p spectral function and a sizable suppression of the peak strength. Because the 
cross sections from the collision rate are evaluated on the p meson mass shell, the imaginary self 
energy of p which comes from the collisions does not vary with M. Therefore the spectral function 
doesn't vanish when the two-pion threshold is reached. 

V. SUMMARY 

Based on the idea that in thermal systems a particle can be captured and can also be produced 
by the thermal background, and the fact that the collision rate is the difference of the capture and 
the production rate, we have calculated the p decay width at finite temperature and the collision 
rates of the p due to the pn scattering and pN scattering in hot and dense hadronic matter. In 
high temperature and/or high density hadronic matter, the collision rate is much larger than the 
decay width correction due to the one-loop self-energy modifications. For the collision rates, the 
contribution from pir collisions is the most important one in the high temperature pion gas, while 
at low temperatures and high density nuclear matter the pN contribution is more important. 

The collision rate of p mesons with pions uses an effective Lagrangian for the aiprr interaction 



which is tuned to hadronic phenomenology 1 21 , 22 1 . Note that in pir — > pir scatterings, the a\ 



intermediate state gives the largest contribution to the cross section and to the collision rate. Our 
result for pN scattering is in qualitative agreement with the width increase found in [ psfl , 
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We have shown that the width of the p is increased drastically through its interaction with 
pions and nucleons present in typical heavy-ion collisions. It will be even more increased if finite 
pion chemical potentials are introduced or if other scattering partners are considered. This result 
will be important not only for the dilepton spectra, but also for the general dynamics of heavy ion 
collisions. We believe that our results adds to the consensus building in this direction. 
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Abstract 

The modification of the width of rho mesons due to in-medium decays and collisions 
is evaluated. The decay width is calculated from the imaginary part of the one-loop self- 
energy at finite temperature. The collision width is related to the cross sections of the p + 
pion and the p + nucleon reactions. A calculation based on an effective Lagrangian shows 
the importance of including the direct pir — > pit scattering which is dominated by the a,\ 
exchange. A large broadening of the spectral function is found, accompanied by a strength 
suppression at the pole. 



PACS : 25.75.Dw, 14.40Cs, 21.65.+f, 13.75.Lb 
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I. INTRODUCTION 



It is generally believed that ultrarelativistic heavy ion collisions produce matter far from the 
ground state and provide conditions of high temperature and/or high density suitable for investi- 
gating modified particle properties in the medium. The vector mesons, e.g. the rho, are particularly 
important components of the matter produced in ultrarelativistic heavy ion collisions. It has been 
speculated that the restoration of chiral symmetry induces a mass shift of the p meson [1]. The 
direct dilepton decay channel, p — > e + e~, appears to be an ideal probe for studying the in-medium 
properties of the p: a possible mass shift, a modified dispersion relation, or a width modification 
[2-8]. Unfortunately, due to the interplay between those effects, it seems to be difficult to iso- 
late the different effects from invariant mass measurement alone. The in-medium mass shift and 
broadening of the width of a particle have been related to the real part of the forward scattering 
amplitude on the constituents of the medium and thus to the corresponding cross sections. These 
results indicate a small mass shift but a large increase of the width both in intermediate and high 
energy heavy ion collisions [5,9,10]. The goal of the present paper is to verify these results using 
slightly different techniques. We will point out the differences as we proceed. 

Recently, dilepton measurements at the CERN/SPS energy have received quite some attention 
[11]. The e + e~ mass spectra (in S + Au collision at E/A = 200 GeV and in Pb + Au collision 
at E/A = 160 GeV preliminary data) seem to indicate a p-peak suppression in central collisions. 
The absence of p-peak may be due to a large width of the p meson as compared with its free value. 
Of course, a confirmation of such a conjecture must include a dynamical simulation for the space 
time evolution of the colliding system. In any case, it is important to quantitatively address this 
issue independently. 

The in-medium effects on the p meson have been widely investigated using finite temperature 
field theory at the one-loop level [3,12,13], but the most important source for the width modification 
occurs at the two- loop level [6]. The present paper focuses on the p meson width modification up 
to the two-loop level. The pole shift of the p, as emerged from many-body calculations is small 
[5,8-10]. 
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The in-medium properties of the p are due to its interaction with the constituents of the 
medium. In the hadronic matter produced in ultrarelativistic heavy ion collisions, pions and 
nucleons are important constituents, so we take into account both pn and p N reactions in our p 
width calculations. Interactions with other hadrons are expected to increase the width above the 
values found here. 

The paper is organized as follows: In section II, we present the general formulae for the collision 
rate of in-medium p's. We present a calculation of the thermal averaged in-medium width. In 
section III, we discuss the cross sections for p tt- and pN- reactions from which the collision rates 
are determined. The width modifications due to pn- and pN- reactions are calculated at finite 
temperature and density in section IV. Section V concludes with a summary. 



II. FORMALISM 

The principal elementary reactions which tend to thermalize p's in hot and dense matter are 
the channels p\ tt 2 — > p 3 vr4 and/or p\ N 2 — > p 3 N4 . A p with arbitrary momentum p\ (and energy 
lo) can be captured by the thermal background. The rate for this direct capture [14] is: 

Td ^ ) = h I dn ME2)(l±nMj)(l±n 4 (E 4 j)\M(l2^M)\ 2 , (2.1) 

A p with momentum p\ can also be produced from the thermal background by the inverse 
reaction: p 3 714 — > p\ r K 2 . This inverse rate (omitted in many similar treatments) is 

= ±-j<Kl{l± n 2 (E 2 )) n 3 (E 3 ) n 4 {E 4 ) \M{M -> 12)| 2 , (2.2) 



where 



and 



dQ, = dp 2 dp 3 dp 4 (2tt) 4 S(p 1 + p 2 - p 3 - Pa) , (2.3) 



dpi 



(2tt) 3 2Et ' 

ni(Ei) are Bose-Einstein or Fermi-Dirac distribution functions, '±' corresponds to Bosons and 
Fermions, respectively. 



The total collision rate for bosons is the difference 



r C0 » = r d M-riM. (2.4) 

This expression can be derived from quantum field theory [15]. From finite temperature quan- 
tum field theory, the width of any particle, T(lo), is given by the imaginary part of its self-energy, 
Imn(cj), 

ImU(io) = -luT(uj) , (2.5) 

which provides a link between field theory and statistical mechanics. 
The ratio of T^lo) to Ti(u) is an universal function of to [16], 



exp(u;/T) , (2.6) 



because the single-particle distribution functions satisfy 

= exp(a;/T) (Boson) , 

n B {io) 

1 -np{io) = exp ^ / T ) (Fermion) . (2.7) 
Thus the total collision rate is given by: 

r co11 = (1 - e-P") T d {oj) . (2.8) 

In the simplest scenario of equilibrated hot and dense matter, we can calculate the thermal 
average of the collision rate: 

^coii Jd 3 pir co11 m (w) J d 3 Pl T d e-P" 



J d 3 pini(cj) / d 3 pini(cj) 
Using (2.1) and (2.9), the thermal rate can be written as 

=coll_M^2 f°°, T w 2 2 



(2.9) 



r C °" = 1 2 / ds-r— ^ r ^A(s,m?,m|)X(v / s,mi,m 2 ,T,// i )cri2-v34(s), (2.10) 



Pl J S0 2(2vr) 4 ^ 

where 



Pl=M /^" lM ^ 



(2tt) 



is the number density of the p's we considered, 0-12^.34(5) is the cross section of pn (or pN) 
collisions, and so = (mi + m2) 2 , Afi are the spin-isospin degeneracy factors. The function K 
assumes different forms for different statistics [17]. Below we apply classical Boltzmann statistics, 
where 

k = K l (^- S /T)^[(p l +p 2 )lT], (2.12) 

with K\ being the modified Bessel function. p\ , p 2 are the chemical potentials of p and tt (or 
nucleon) respectively. Prom (2.10) the collision rate of p in the hot and dense matter can be 
obtained if the cross sections of the p tt- or pN- processes are known. 



III. p?T AND pN CROSS SECTIONS 

For p tt collisions, we use both the scattering amplitude of the elementary processes method 
and the resonance contribution approximation to calculate the cross sections a p7T . We consider the 
scattering amplitudes for the elementary reactions in Fig. 1 which can be obtained by cutting the 
two- loop diagrams for the p meson self-energy as demonstrated in Ref. [18]. 

To determine the p N cross section <j p n, a method based on the resonance model is used to get 
a lower estimate. 



A. pn cross section 



The standard Breit-Wigner formula for p tt cross section is 



t x ivs- m R f + ry 4 

Here 'R' refers to the intermediate resonance, and Br is the branching ratio of its decay into tt p. 
V^s is the total center-of-mass energy, q is the three-momentum of the p meson in the resonance 
frame, 

q = -L\V 2 ( S: ml,m 2 p ) : (3.2) 
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with the kinematical triangle function X(x,y,z) = x 2 — 2x(y + z) + (y — z) 2 . F s and Fi are spin 



and isospin degeneracy factors: 

2Ir + 1 



2Jr + 1 



(3-3) 



(2I n + l)(2/ p + l) ' (2 S7r + l)(2 Sp + l) ' 

We take into account the following resonances: 0(1020), ai(1260), a 2 (1320), a/(1420), vr(1670) 
with the properties listed in [19]. 




(2) 



(5) 





(2) 



FIG. 1. The Feynman diagrams, (a) for p ^ —> p ^ reactions and (b) for p°tt° — > p°-K° 
reactions. 

The pion scattering amplitude on any hadronic target vanishes at zero pion energy in the target 
rest frame in the zero pion mass limit (Adler's theorem). Following Ref. [9] one can include an 
additional factor to the resonance contributions to a pn , 
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which is normalized to 1 at s = m 2 R ; for s > m R this factor is taken to be one. The resulting 
resonant cross section a pv is shown in Fig. 2(a) (dashed curve). 

The cross section for pir scattering can be calculated by using an effective Lagrangian for the 
hadronic interaction pir — > pir. The scattering proceeds through s- and t- channels respectively, 
where the 7r, p, lu, </>, a\ and u/(1420) might be intermediate states. Those processes are constructed 
from 7T7T/9, 7rpcj and npai vertices. The effective interaction Lagrangians we use are the following 
[20]: 

tp-K-K = igp-Knin x d^n) ■ p^ , (3.5a) 

£ PP ™ = g%A^ x Pn) ■ x p") , ( 3 - 5b ) 

£pvuj = gpTToj^^d^UJ^daPr ■ if . (3.5c) 



The La grangian Cpna\ is too long to be presented here [21,22] but for the physical a\(h) — )■ 
p(q) 7T (p) decay, for example, the vertex function is 

fV = i{fAg/iu + gAluK + hAq^q,,) , (3.6) 

where 

Ja = -j= [- W + (m - m)k ■ q] , 
gA = _? ?2), 

h A = -^=r)i . (3.7) 

The third term in (3.6) actually does not contribute to the decay width. We chose parameters 
which reproduce the measured a\ — > pir decay width and D/S ratio (g = 6.45, iji = 2.39 GeV -1 , 
r]2 = 1.94 GeV -1 ). The coupling constant g pnv = 6.05 is determined from a fit to the free p — > tttt 
decay width, g p7ru = 12.40 GeV -1 is determined by the cj — > pir — > 7T7 decay [20]. The coupling 



constant g p7R y = 3.55 GeV -1 is determined by the lo' — > pn decay. Note that we believe that it is 
important for our effective Lagrangians to generate good phenomenology at energy scales relevant 
for the applications considered here. 

The Feynman diagrams for p tt scattering are given in Fig. 1. The cross section is then calculated 



as 



where 



t ± = ml + m 2 p - ^-[(s + m%- m 2 p ) (s - m 2 n + mjj)] ± ^A(a, ml,m 2 p ) . (3.9) 

For p° tt^ —> p° tt^ scattering, the squared amplitude is obtained from the processes shown in 
Fig. 1(a) with the appropriate interference contribution: 

1 



\ M \ =7^ iT^ -^1^1+^2+^3+^4+^51% (3-10) 

(2s v + lj(2s p + 1) 

which is sum over the spins of the initial states as well as over those of the final states. The 
p° tt^ —> p^ 7T° scattering is also considered, The squared amplitude has the same form as (3.10) 



but no contributions from A^2 and M.§. For p° n° scattering, Fig. 1(b), \M.\ can be calculated 
similarly. 

Due to the substructure of hadrons, for vertices in the t-channel Feynman diagrams we use 
form factors 

^ = ^T> ( 3 - n ) 

where a indicates the species of the exchange particle. A = 1.8 GeV is taken for all vertices. 
The numerical integration for the reaction pn — > n — > pn (t-channel) results in a singularity. We 
regulate it with the effective approach of Peierls [23]. This technique offers a pragmatic solution to 
the possible unitarity violation at high energies that can result from the use of tree-level diagrams. 
A more complete solution would perhaps comprise a K-matrix multichannel calculation. However, 
such an endeavor at finite temperature and finite density also carries ambiguities of its own. 
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After computing the amplitude, the 7rp scattering cross section is obtained from (3.8). The 
result is shown in Fig. 2(a) (solid curve). The resulting pir cross sections indicate that the difference 
between the resonance model and the scattering amplitude method becomes large at high energies. 
The resonance model may be valid only at low energies. Our effective Lagrangian technique also 
treats below-threshold resonances and also receives a contribution from t-channel diagrams and all 
their associated quantum interferences. 




FIG. 2. The cross sections, (a) for pn — > pn reactions, and (b) for p N — > p N reactions. 



B. pN cross section 

The pN cross section can be calculated within the resonance model as [9,10], 

Here 

q =^L\^( aj rn%,m 2 p ), (3.13) 

and = trip + mp?. The summation is performed over all baryonic resonances with masses 
above the pN threshold and below 2200 MeV as quoted in [19]. The result of the pN cross section 
calculated from (3.12) is shown in Fig. 2(b) (dashed curve). 



The p N cross section obtained within the resonance model may be valid only for low energies, 
while at high energies one should calculate it from the quark model. Furthermore, it has been 
shown recently that the resonances below the p N threshold, like the N*(1520) play an important 
role in the p-nucleon dynamics [24], so in the baryonic sector we adopt the pN cross section as 
determined by Kondratyuk et al. Ref [25], which includes also resonances below the pN threshold. 
The numerical result is shown in Fig. 2(b) (solid curve). 

IV. BROADENING THE WIDTH OF THE RHO AT FINITE TEMPERATURES 

AND DENSITIES 

The modified properties of hadrons in hadronic matter, e.g. the width increase of the particles, 
are due to their interactions with the constituents of the medium. We assume that ultrarelativistic 
heavy ion collisions generate a hot and dense hadron gas of pions and nucleons and consider the 
p mesons interacting with this thermal gas. We realize that this is a simplification of the complex 
dynamics which will be generated in all collision of heavy ions. Our results should be taken as 
lower estimates of the realistic effects. 

As discussed in section II, the rho meson width is related to the imaginary part of the rho 
self-energy. The total width of the p meson is the physical decay width plus the collision width, 

r total = pdecay + ^coll. _ ^ 

The decay width is related to the imaginary part of the rho self-energy at the one-loop level. 
The collision width is related to the imaginary part of the rho self-energy at the two-loop level. 
Let us first calculate the decay width of the p mesons at finite temperature. 

A. The decay width of the rho at finite temperature 

To evaluate the p meson self-energy from the tttt loop we consider the pirn interaction La- 
grangian from the previous section. Using the Feynman rules of thermal field dynamics, the 
temperature dependent rho meson self-energy is [12] 
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lWg,T) 



9piriT 



J d*k {-2 9ti J(k 2 - ml)n n (k) + (q + 2k)^q + 2k) v - 



(2vr) 3 

5(k 2 — ml) 



,,s S((k + q) 2 -ml) . 
-n*{k) + " „ n v {k + q) 



(4.2) 



(A; + q) 2 — m 2 k 2 — m 2 

Because of the lack of Lorentz invariance in the medium [3,12,26], the p self-energy at finite 
temperature separates into a transverse part, Ht(q, T), and a longitudinal part, Hl(q, T). In the 
rho meson rest frame, 



U(co,q^0,T) = U L (u,q^0,T) = U T (u, q ^ 0,T) 



(4.3) 



thus the real part of the self-energy is, 



ReU(co,q^O,T) = - 9 ^-jdk^- ^3 



4k 2 



4E 2 



n n (E k ) 



(4.4) 



where E k = \Jk 2 + m 2 . The real self-energy will lead to a mass shift of the p at finite temperature, 
which is determined by the solution of 



u? -m 2 p + ReU(uj,q^0,T) =0. 

It is known from previous calculations that this mass shift is small [3]. 

The imaginary part of the self-energy can be obtained by using cutting rules [16], 
a 2 f h 4 

lmU(co,T) = / dk (2n n (E k ) + 1) [6{lj + 2E k ) - 6{lj - 2E k )} . 

D7T J £j k 

From ImII(cj,T), the corresponding width of the p — > tttt decay follows to be 

r deca y(T) = Al_ {u) 2 _ im 2 )S/ 2 



(4.5) 



(4.6) 



2n,(|) + l 



(4.7) 



This decay width varies with the temperature as shown in Fig. 3. Obviously the decay width 
increases with temperature as well as with the pion chemical potential. For T = 150 MeV, when 
^ = 0, Ar p = r* eca y(T) - r^ ca y(0) = 25 MeV, while for ^ = 135 MeV, Ar p = 70 MeV. 
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FIG. 3. The p — > tttt decay width as 
a function of temperature T at two fixed 
values of the pion chemical potential. 
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FIG. 4. The p collision rate vs. the tem- 
perature due to p 7r interaction in two dif- 
ferent approaches. 



B. Collision rate of rho at finite temperature and density 

From (2.10) and with the cross sections computed in the previous section, the collision rates of 
the p mesons can be calculated numerically. 

First the rho collision rate is calculated due to pn reactions in the finite temperature meson 
medium. The contributions from p ^ —> p ^ , p ^ — > p^ir and p ir° — > p ir° processes are 
taken into account. Fig. 4 shows the the collision rates as calculated from the scattering amplitude 
calculation and from the Breit- Wigner resonance approximation at two values of the pion chemical 
potential, p, n = and p, v = 135 MeV. For the same temperature and chemical potentials and using 
the p-K cross section from the scattering amplitude method, one can get a larger collision rate of p 
mesons. 

The collision rate from p N reactions is plotted in Fig. 5. In Fig. 5(a), we use T = 100, 150 
MeV and show the rate as a function of the nucleon density ratio. In Fig. 5(b), we show the rate 
as a function of the temperature for two values of the nucleon chemical potential, Hn = 0.4, 0.6 
GeV. From Fig. 5(a), we can see the collision rate increasing with the nucleon density for a given 
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temperature, but decreasing with temperature for a fixed nucleon density. The latter decrease is 
due to the fact that the nucleon chemical potential decreases with temperature at a fixed nucleon 
density. Again we find that for given temperature and chemical potential the pure resonance 
contribution yields a lower collision rate than the resonance + quark model. 
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FIG. 5. The collision rate of the p due to p N interactions, vs. (a) the nucleon density for vary 
temperature, (b) the temperature for vary baryon chemical potential. 

We employ the resonance model modified to include the resonances below the pN threshold, 
such as N*(1520), as described previously. The total pN cross section can be evaluated as a function 
of the invariant collision energy yfs and the invariant mass of the p meson M p . In this case, the 
collision width has to be integrated over the spectral function of the p meson as 



r 



p^ 



M X M 2 



ds- 



ifs—mi* 



dM p M p A(M p ) \(s, M 2 m 2 N ) K a pV (s, M p ) , (4. 



Pl Js > 2(2^) 4 ^ j 2m . 
here s' Q = (2m v + m^) 2 . A(M p ) is the spectral function of the p meson in free space taken as 

1 m p T p (M p 



A(M p 



(4.9) 



p> it {Mj - mj) 2 + mjTj ' 

where m p = 770 MeV, T p (M p ) is the mass-dependent width of the p meson [25]. For the same 
process of pN scattering, the collision width of the p meson in this calculation is about 10% larger 
than that from the on mass-shell p meson calculations. 
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In summary, our results indicate that at high temperature and low nucleon density, the collision 
rate is dominated by p + meson reactions, while at low temperature and high nucleon density it 
is dominated by p + nucleon reactions. For example, when T = 100 MeV, p v = 0, pjy = Po, 
r p7r coU - = 12 MeV, while T pN c ° 11 = 180MeV. 

For conditions typical of high energy heavy- ion collisions [5,27], T = 150 MeV, p v = and 
Pn = 0.4 GeV, (pn = 0.39po ) the total width of rho mesons is 

ptotal _ pdecay _|_ p co11 - _|_ p ^ co11 - 

sa 176 + 58 + 63 = 297MeV . 
This dramatic increase of the p width is most certainly important. 




FIG. 6. The spectral function of the p vs. the invariant mass M of the p for fixed values of 
momentum as indicated, (a) In pure pion gas. (b) In the matter of pions and nucleons. 

The spectral function of rho meson is given by 

2ImII 

P(M) = ~ (a; 2 - q 2 - mj - Ren) 2 + (IuiII) 2 ' (4 ' 10) 
where q is the three-momentum of the p meson. The vacuum part of the self-energy II, which is 



given in [28], can only depend on the invariant mass M = Vw 2 - q 2 . We plot the p meson spectral 
function in fig. 6(a) for a pure pion gas at a temperature T = 150 MeV and in fig. 6(b) for the 
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matter of pions and nucleons at T = 100 MeV and nucleon density pn = Po = 0.16 fm -3 . We find 
a broadening of the p spectral function and a sizable suppression of the peak strength. Because the 
cross sections from the collision rate are evaluated on the p meson mass shell, the imaginary self 
energy of p which comes from the collisions does not vary with M. Therefore the spectral function 
doesn't vanish when the two-pion threshold is reached. 

V. SUMMARY 

Based on the idea that in thermal systems a particle can be captured and can also be produced 
by the thermal background, and the fact that the collision rate is the difference of the capture and 
the production rate, we have calculated the p decay width at finite temperature and the collision 
rates of the p due to the p tt scattering and pN scattering in hot and dense hadronic matter. In 
high temperature and/or high density hadronic matter, the collision rate is much larger than the 
decay width correction due to the one-loop self-energy modifications. For the collision rates, the 
contribution from p n collisions is the most important one in the high temperature pion gas, while 
at low temperatures and high density nuclear matter the pN contribution is more important. 

The collision rate of p mesons with pions uses an effective Lagrangian for the aipn interaction 
which is tuned to hadronic phenomenology [21,22]. Note that in pn — > pn scatterings, the a\ 
intermediate state gives the largest contribution to the cross section and to the collision rate. Our 
result for pN scattering is in qualitative agreement with the width increase found in [28]. 

We have shown that the width of the p is increased drastically through its interaction with 
pions and nucleons present in typical heavy-ion collisions. It will be even more increased if finite 
pion chemical potentials are introduced or if other scattering partners are considered. This result 
will be important not only for the dilepton spectra, but also for the general dynamics of heavy ion 
collisions. We believe that our results adds to the consensus building in this direction. 
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